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This article describes a statistical model to analyze the signal intensity received at 
the solid-state imaging (SSI) camera of the Galileo optical communications system 
from an Earth-based transmitter (GOPEX) demonstration. The analytical model 
assumes that the optical beam possesses a Gaussian profile and the communication 
channel has a log-normal scattering characteristic. The atmospheric-induced jitter is 
modelled as two independent zero mean Gaussian random variables. By modelling 
the system parameters as a set of independent and identically distributed (iid) 
random variables, the combined impact of uncertainties due to system parameters 
and the turbulent atmosphere can be approximated by a log-normal distributed 
signal intensity at the spacecraft. A Monte-Carlo software simulation package has 
also been developed to compute the confidence interval probabilities for general 
optical beam profiles. Numerical results show that the approximation is valid for a 
wide range of operation scenarios. 


I. Introduction 

In this article, analytical expressions for the probabil- 
ity density functions (pdf’s) of the optical field intensity 
and the number of observed photoelectrons per pixel of 
Galileo’s solid-state camera for a laser pulse interval in the 
Galileo optical communication from an Earth-based trans- 
mitter (GOPEX) experiment are derived. These pdf’s play 
a critical role in assessing the performance of this optical 
uplink and provide a means of acquiring meaningful design 
predicts. Furthermore, the aforementioned pdf’s shed light 
on the characteristics of the underlying processes involved 


in this experiment and the impacts of various parameters 
on the detection probability. 

It is assumed that the main sources of disturbance 
are the atmosphere-induced jitter, for the short-exposure 
model at hand, and the log-normal large-particle channel 
scattering, which proves to be significant for the operating 
characteristics of this experiment. In general, there are a 
number of parameters, described below, that directly im- 
pact the observed signal intensity at the spacecraft. These 
parameters are assumed to be independent and identically 


distributed (iid) random variables. It is imperative to note 
the importance of this observation, since as a result of 
this assumption the signal intensity in the absence of at- 
mospheric effects may be approximated by a log-normal 
random variable as well. In this article, the log-normal 
channel scattering is shown to be the dominant effect for 
the operating characteristics of this experiment. Thus, 
it is concluded that the combined impact of random sys- 
tem parameters and turbulent atmosphere would yield a 
log-normal distributed signal intensity at the spacecraft. 
Before going any further, a description of the system pa- 
rameters that directly impact the observed signal intensity 
and which are needed to conduct the ensuing analytical 
study of the GOPEX uplink is presented. 

The optical field intensity at the spacecraft is a func- 
tion of transmitter optics efficiency, ty , receiver optics ef- 
ficiency, rj r , atmospheric transmittance, T) atm , transmitter 
laser energy (in joules), W t , distance from the spacecraft, 
Z, and the predefined angular beam diameter, 0 S . Since 
the number of photoelectrons per pixel, observed over a 
laser pulse duration, is of importance here, one must also 
include the effective detector area (pixel area), A rec , and 
the quantum efficiency of the photodetector, r/ c , in deter- 
mining the overall statistical modelling. In the event of a 
uniform beam pattern, and for Z >> 1, one can describe 
the optical field intensity in joules/m 2 at the spacecraft as 


optical field intensity at (x,y) in the plane of observation 
may be approximated by 
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where now I t accounts for the log-normal atmospheric 
scattering effect, r]o is the atmospheric loss at zenith, 
and 9 is the spacecraft’s zenith angle. This implies that 
rjatm = , which agrees with most experimental and 

theoretical studies of atmospheric absorption. In this anal- 
ysis, assume a negligible pointing error, and thus the field 
intensity must be evaluated for (x,y) = (0,0). Further- 
more, since the receiving area of the charged-coupled de- 
vice (CCD) camera is significantly smaller than the beam 
footprint at the spacecraft location, the CCD camera may 
be considered a point detector. Therefore, Eq. (3) at 
(z,2/) = (0,0) accurately describes the field intensity for 
all the detectors of the CCD camera. 

The pdf of 7/ is given by the well-known log-normal 
density 
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where < I(x,y,t) > t describes the intensity in the plane 
of observation at coordinates (x,y) } and < . > ( signifies 
short-term intensity averaging over time, since the obser- 
vation interval is limited to only 6 to 12 nsec, the laser 
pulse duration. The number of observed photoelectrons in 
the laser pulse duration is then 
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where crj is, in turn, given by [1] 
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These expressions will be modified in the next section to in- 
clude the Gaussian beam profile and the log-normal large- 
particle scattering effect. 


II. Uplink Analysis: Gaussian Beam Profile 

In this section, a transmitted beam with a Gaussian 
profile is assumed. Therefore, for a short exposure model, 
one may assume that the Gaussian beam profile maintains 
its integrity in the plane of observation. In this case, the 
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In the above equation, A is the wavelength of the laser 
in meters, H is the height of the atmosphere, and C n (h) 
is the medium index of the refraction structure constant. 
Major stumbling blocks in determining an accurate esti- 
mate of af are the dependency of C n on various random 
channel parameters and the unavailability of an accurate 
model for C n . However, in the literature a number of ap- 
proximate models for this structure constant are available. 
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For a complete list, refer to [2]. Among those references, 
the following expression, originally proposed by Hufnagel, 
approximates this structure constant: 

C 2 (h) = 8.2 x 10“ 56 V r2 /i 10 e“ ;i/1000 

-f2.7 x io~ 16 e" h / 1500 (6) 


where k is the wave number. For the parameters of this 
experiment, a is estimated to be 9 ^rad for near zenith and 
14 ^rad for a zenith angle of 60 deg. Define the random 
variable I g as follows: 
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where h is the altitude in meters above sea level and V is 
the rms wind velocity, averaged over 5 to 20 km altitude, 
and is considered to be Gaussian distributed about a mean 
of 27 m/sec with a standard deviation of 9 m/sec. This 
model, however, is valid only for altitudes in excess of 5 
km. For altitudes below 5 km, an approximate expression 
from Hufnagel and Stanely for the structure constant is 
given by 


C 2 (h) » 1.5 x I0~ l3 /h ( 7 ) 

For the parameters of the GOPEX project, and based 
upon the above approximations, <t\ < 0.9. 

In a turbulent medium, a Gaussian beam, as described 
above, experiences deflection in various directions due 
to gaseous blobs and other turbulent particles flowing 
through the path of the beam. For a short exposure model 
(6- to 12-nsec observation interval), the atmospherically 
induced deflections may accurately be modelled as atmo- 
spherically induced pointing “jitter.” This implies that x 
and y in Eq. (3) may now be viewed as two independent 
and zero mean Gaussian random variables with standard 
deviation <rZ. Thus, a represents the rms angular (half- 
beam) jitter due to turbulence. When pointing error is 
present, its standard deviation may be added directly to 
a. However, caution must be exercised in applying the 
ensuing results to a model with a constant pointing error. 
In this event, x and y are nonzero mean Gaussian random 
variables. For this analysis, however, the constant point- 
ing error is considered to be negligible, and thus one may 
consider x and y as zero mean Gaussian random variables. 

The standard deviation a can be estimated for a given 
index of refraction profile by using the following expression 
[3]: 
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It can readily be shown that for jointly Gaussian x and 
y, Ig is a special case of beta distributed random variables 
with pdf 
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where /? = 0 2 /8<r 2 . In the following, (3 ranges from 5 to 19 
(6 S = 110/irad and a « 9 to 14/zrad) for all practical pur- 
poses. It is interesting to note that for (3 — 1 this density 
reduces to a uniform density, signifying the detrimental 
impact of a turbulent medium. 

The pdf of It — Igh can now be found. It represents 
the combined impact of log-normal medium scattering and 
atmospherically induced pointing jitter. This pdf is ex- 
pressed as 
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where Q(x) = / r °° exp (s 2 /2) ds. For a reason- 

able range of system parameters, this density is depicted 
in Figs. 1 through 4. It is important to note that ai plays 
a critical role in determining the behavior of this random 
variable. Unfortunately, for zenith angles in the range of 
40 to 55 deg, A = 0.532 ^m, and a rms wind velocity of 
27 m/s, <r; proves to be in excess of 0.2, resulting in a pdf 
that can accurately be approximated by a log-normal den- 
sity function. This implies that for all practical purposes, 
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These expressions can be numerically evaluated for a de- 
sired set of system parameters. The validity of approxi- 
mating the characteristics of I t with that of a log-normal 
random variable is examined in a number of plots (see 
Figs. 5 through 8 ). As noted, for a wide range of system 
parameters the penalty for this approximation is comfort- 
ably small. 

One can now express the number of observed photo- 
electrons in a pulse interval as 
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where r, iot = flLi Vi, with m = Vr, K = Vt, Vs = vT W , 
= Vc, and 775 = A r6C . If one considers 77 , ’s as iid random 
variables, it is possible to approximate r) tot with a log- 
normal random variable. Due to a lack of sufficient data 
for accurate characterization of these parameters or for 
considering the worst-case scenario, one may assume that 
77 ,- is uniform over ( 77 ™ 
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where m„ = £? =1 m„ t and <r 2 = £? =1 <r 2 . with 
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Because the product of two log-normal random vari- 
ables yields log-normal statistics 
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where N toi = Vtothy wn = m v - h rn t , and a 2 N = cr 2 -j- er 2 . 
Finally, the total number of photoelectrons observed over 
a pulse interval is given by 
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Note that an alternate means of computing <j 1 v and m/v 
is to first compute the mean and mean-square of r} tot I ti 
which are given by 
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and use the following expressions: 
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<r% = log [E{N? ot }/E 2 {N tot } + 1] (23) 

mjv = log [E{N tot }] — -<t 2 n (24) 

It is imperative to note that the accuracy of the above 
model is highly sensitive to cr\ and other system parame- 
ters. In particular, caution must be exercised in assuming 
a log-normal scattering model for values of crj > 0.75, since 
for such values of < 7 / the variance of the normalized in- 
tensity fluctuation due to scattering (i.e., e a ' - 1) exceeds 
0.75. In this event, log-normal statistics are no longer valid 
in describing the characteristics of the scattering channel. 
Instead, an exponential pdf must be employed to char- 
acterize the scattering channel [3,4]. For the experiment 
at hand, < 7 / for V < 27 m/sec and 9 < 55 deg is less 
than 0.9. This, therefore, makes the log-normal assump- 
tion rather suspect for some portion of this experiment 
where the spacecraft takes on large zenith angles. How- 
ever, since the rms wind velocity is a random variable, and 

is highly sensitive to this parameter, one may approx- 
imate the channel scattering effect with the log-normal 
statistic. 

III. Uplink Analysis: General Beam Profile 

This section provides a brief description of a software 
simulation package that was developed for the analysis of 
the GOPEX uplink when the optical beam possesses a gen- 
eral profile. To compute confidence interval probabilities, 
one has to resort to a Monte Carlo simulation to include 
the impacts of atmospherically induced jitter, log-normal 
channel scattering, and uncertainties in other system pa- 
rameters that are modelled as uniformly distributed ran- 
dom variables (see above). 

This package reads the beam matrix profile, recorded 
at a known distance, and evaluates the beam profile at 
the location of the spacecraft. The turbulent medium is, 
once again, modelled as a zero mean Gaussian pointing 
jitter in both x and y coordinates. This program provides 
the mean, standard deviation, and the pdf of the num- 
ber of observed photoelectrons in a pulse duration for a 
specified interval. The confidence probability for a given 
confidence interval is also computed. Since no theoreti- 
cal results were available to test the validity of the results 


of this program for a general beam profile, a Gaussian 
beam was specified as the input beam profile. The system 
parameters were set at the following values: W t = 0.25 
joules, A = 0.532 //m, 9$ = 110/irad, rjt — 0.705, r) r = 0.24, 
tj c = 0.37, rjo = 0.715 (atmospheric attenuation at zenith 
and station altitude of 2.286 km [5]), A rec = 0.01825 m , 
Z - 0.6 x 10 6 km, V = 27 m/sec, and 9 = 55 deg. From 
the above analytical results for the Gaussian beam profile, 
the mean and standard deviation of the observed photo- 
electrons were estimated at 111.991 x 10 3 and 126.08 x 10 3 , 
respectively. The large standard deviation is a clear indi- 
cation of log-normal statistics due to dominant channel 
scattering. The simulation for 10,000 samples yielded the 
following mean and standard deviation: 110.834 x 10 3 and 
119.535 x 10 3 . These numbers improved to 111.950 x 10 3 
and 124.361 x 10 3 for 50,000 samples. It is quite clear 
that a reasonable accuracy may be achieved with a rela- 
tively small number of samples. However, with the aid of 
stratified or importance sampling, the required number of 
samples may be drastically reduced. The above calculation 
was repeated for a 60-deg elevation, i.e., 9 = 30 deg, with 
the remaining parameters fixed at the above values. From 
theoretical results, the mean and standard deviation were 
found to be 142.025 x 10 3 and 97.166 x 10 3 , respectively. 
For 10,000 samples, simulation yielded 141.27 x 10 3 and 
94.877 x 10 3 as the mean and standard deviation of the 
observed photoelectrons, respectively. Once again, these 
numbers improved to 142.095 x 10 3 and 96.536 x 10 for 
50,000 samples, clearly indicating the consistency of the 
simulation. 


IV. Conclusion 

This article described a statistical model to evaluate the 
signal intensity received at the solid-state imaging camera 
for the GOPEX demonstration. The model includes the 
efTect of log-normal channel scattering, atmospherically 
induced jitter, and uncertain system parameters. It has 
been shown that the resulting probability density function 
can be well approximated by a log-normal distribution. A 
Monte Carlo simulation software package has been devel- 
oped to analyze the uplink performance when the opti- 
cal beam profile is non-Gaussian. The confidence-interval 
probability can be computed to analyze the GOPEX ex- 
perimental data during demonstration. 
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PROBABILITY DENSITY FUNCTION 



Fig 1. Probability density function of the normalized intensity, 
/3 = 3. 



Fig. 3. Probability density function of the normalized Intensity, 
(3 = 7. 



Fig. 2. Probability density function of the normalized intensity, 
(3 = 4.6. 



Fig. 4. Probability density function of the normalized intensity, 
(3 = 10 . 
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PROBABILITY DENSITY FUNCTION 



Fig. 5. Probability density function of the normalized intensity, 
P = 5 and a i = 0.5, and its approximated log-normal counter- 
part. 



Fig. 7. Probability density function of the normalized intensity, 
P = 5 and <jl = 0.2, and its approximated log-normal counter- 
part. 
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NORMALIZED INTENSITY 

Fig. 6. Probability density function of the normalized Intensity, 

P = 10 and <j i = 0.2, and its approximated log-normal counter- 
part. 



Fig. 8. Probability density function of the normalized intensity, 
P = 5 and cr^ = 0.5, and its approximated log-normal counter- 
part. 
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